Collective Light Emission of a Finite Size Atomic Chain 
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Radiative properties of collective electronic states in a one dimensional atomic chain are investi- 
gated. Radiative corrections are included with emphasize put on the effect of the chain size through 
the dependence on both the number of atoms and the lattice constant. The damping rates of collec- 
tive states are calculated in considering radiative effects for different values of the lattice constant 
relative to the atomic transition wave length. Especially the symmetric state damping rate as a 
function of the number of the atoms is derived. The emission pattern off a finite linear chain is also 
presented. The results can be adopted for any chain of active material, e.g., a chain of semiconductor 
quantum dots or organic molecules on a linear matrix. 

PACS numbers: 37.10.Jk, 42.50.-p, 71.35.-y 



I. INTRODUCTION 

Optical lattice ultracold atoms continue to be of in- 
terest for more and more researches of different branches 
of physics [l|. Big attention is given for the realization 
of different condensed matter models that provide a test 
system for achieving a deep understanding of fundamen- 
tal physics and answering open questions in the subject 
[2|, beside their applications for quantum information 
processing Q. In general, the main objective is to con- 
sider optical lattice ultracold atoms as artificial crystals 
with a wide range of controllable parameters. 

Optical lattices form of counter propagating laser 
beams to get standing waves in which ground state ultra- 
cold atoms are loaded 4, 5]. The atoms experience opti- 
cal lattice potential with lattice constant of half wave 
length of the laser. Low dimensional lattices can be 
achieved with different geometric structures and symme- 
tries ja]. In conventional solid crystals the lattice con- 
stant and the symmetry of the lattice is fixed through 
the different chemical bonds that responsible for the for- 
mation of the crystal. The advantage of optical lattices 
is due to the controllability of the lattice constant and 
symmetry through controlling the external laser field [l| . 

Collective states of electronic excitations play a central 
rule in solid crystals and molecular clusters and they usu- 
ally termed excitons 0, [8j . They induced by electrostatic 
interactions among the lattice atoms or molecules, where 
an electronic excitation can be delocalized in the crystal 
through energy transfer. Collective states can dominate 
the electrical and optical properties of the material, and 
especially they strongly affect the excitation lifetimes and 
give rise to dark and superradiant states. In such mate- 
rial the lattice constant is few angstroms which is much 
smaller than the electronic transition wavelength, and 
hence one can use electrostatic interactions, e.g. reso- 
nance dipole-dipole interactions, and to neglect radiative 
corrections altogether. 

Electronic excitations in optical lattice ultracold atoms 
are of big importance, e.g., for optical lattice clocks [9j, 
and for optical lattice Rydberg atoms [ly]. In our previ- 
ous work we introduced excitons for optical lattice ultra- 



cold atoms in one and two dimensional set-ups [111 . Il2| . 
We concentrated mainly in the Mott insulator phase with 
one and two atoms per lattice site. We treated both large 
and finite atomic chains [13 - ll5j . and we calculated the 
damping rate of excitons into free space and their emis- 
sion pattern [loHlq . In all of our previous researches we 
exploited electrostatic interactions for the formation of 
collective states, mainly resonance dipole-dipole interac- 
tions. But for typical optical lattices the lattice constant 
is few thousands of angstroms, which can be of the order 
of the electronic transition wavelength, and hence radia- 
tive corrections can be significant. 



In the present paper we investigate a one dimensional 
finite chain of atoms where the lattice constant can take 
any value relative to the atomic transition wavelength. 
Finite atomic chains have been realized recently in a 
number of optical lattice experiments [19|, [2(| . We em- 
phasize the influence of radiative corrections on the for- 
mation of collective sates and their damping rates, where 
we exploit general collective states with emphasize on the 
most symmetric one. We derive the condition for the 
validity of applying electrostatic interactions, which we 
used in our previous work. Few studies treated the collec- 
tive effect on the optical properties of finite atomic chain 
of several atoms [2l|, but extensive study done for two 
atoms in the radiative regime [22j | , and in which we com- 
pare our results. We extract how the damping rate de- 
pends on the chain size, namely on the number of atoms 
in the lattice. Furthermore, we calculate the emission 
pattern off a finite atomic chain. 



The paper is organized as follows: in section 2 we 
present a finite one dimensional atomic chain and dis- 
cuss the energy transfer parameter due to dipole-dipole 
interactions in the radiative regime. Then in section 3 we 
calculate the damping rates for different collective states 
and several chain sizes. The emission pattern for col- 
lective states is calculated in section 4. The summary 
appears in section 5. 



II. FINITE ONE-DIMENSIONAL ATOMIC 
CHAIN 

We consider a finite one dimensional atomic lattice, 
where the number of atoms is N with lattice constant 
a, as seen in figure (1). The atoms are considered to 
be two-level systems with electronic transition energy 
Ea = hujA. An electronic excitation can delocalize in 
the lattice by transferring among the atoms. The elec- 
tronic excitation Hamiltonian is given by 



H„ 



J2 ?iUA B\B n 



y^ hj„ 



iBlBn 



(1) 



where B\ and B n are the creation and annihilation oper- 
ators of an electronic excitation at atom n. For a single 
excitation the operators can be assumed to obey boson 
commutation relations. 




FIG. 1: A finite lattice of TV atoms. The lattice constant is a, 
and the transition dipole /x makes an angle tp with the lattice 
direction. 

The energy transfer among two atoms, n and m, is a 
function of the interatomic distance and given by |23| 



J {lARnm) 




sin(q A R„ m ) COS (tfA-Rnm) 



(QARn 



{qARnmY 



(1-3 cos (p] 

COS jqARnm) 
qAtinm 



(l - COS 2 If) 



(2) 



where the distance between the two atoms is R nm = 
\n — m\a, and /i is the magnitude of the electronic excita- 
tion transition dipole, which makes an angle tp with the 
lattice direction, see figure (1). qA is the atomic transi- 
tion wave number given by Ea — frcqA- Here Ta is the 
single excited atom damping rate 



r A = 



STTEofrc 3 



(3) 



In the limit of A^ > a, where Xa is the atomic transition 
wave length defined by Ea — Hc/Xa-, we can consider 
only energy transfer among nearest neighbor atoms with 
J {qAo) where we take R nm = Q- 

In figure (2) we plot J (qAQ-) I^A as a function of 
qAd for two different polarization directions. Note that 
for typical optical lattice we have Ea — 1 eV, with 
X A « 12405 A, and q A « 4 x 1CT 4 A~ l . For a = 1000 A 
we get qAd ~ 0.5, and ci/Xa ~ 0.08. For ip — 0° we 
obtain J(0.5)/r A w -13.4, and for <p = 90° we get 
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FIG. 2: The scaled interaction J (qAa) /Ta vs. qAd. 
line is for ip = 0° , and the dashed line for tp = 90° . 



The full 



J(0.5)/Ta ~ 5.4. For large qAa the coupling tend to zero 
with oscillations, and the atoms are almost independent. 
In the limit Xa 3> a, or g^a ^C 1, we can ne glec t the ra- 
diative terms (as we did in our previous works |lll4l8l |). to 

get the electrostatic resonance dipolc-dipole interaction 



Jw-- Fa 



4 (gAa) : 



(l — 3 cos 2 ip) 



(4) 



Using the previous numbers, <p = 0° yields J(0.5)/r j 4 w 
-12, and (p = 90° yields J(0.5)/T A « 6, which are 
slightly different from the above exact results. For 
smaller g^a we get much better agreement. In figure 
(3) we plot equations ([2]) and ((4]) for ip = 0°, and in 
figure (4) for tp = 90°. The results justify the use of 
electrostatic dipole-dipole interactions for optical lattice 
ultracold atoms when g^a < 1. 




FIG. 3: The scaled interaction J (qAa) /T a vs. qAa for tp = 
0°. The full line is for equation IJ2)|, and the dashed line for 
equation Q. 
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FIG. 4: The scaled interaction J (g^ia) /I\a vs. g^a for <p — 
90°. The full line is for equation J2]), and the dashed line for 
equation (Q. 



III. COLLECTIVE EXCITATION DAMPING 
RATE 

We start in presenting the free space radiation field 
and its coupling to a finite atomic chain. The free space 
radiation field Hamiltonian is 



H ra d = 2_^Eph(q) « q A a q^' 



(5) 



qA 



where a x and a q A are the creation and annihilation op- 
erators of a photon with wave vector q and polarization 
A, respectively. The photon energy is E p h{q) = hcq. The 
electric field operator is 



qA 



hcq 
^—^ ja q A e qA e- ' - a qA e qA e 



o»q r _ „t o*. 0-iq-T 



}■ 



(6) 

where e qA is the photon polarization unit vector, and V 
is the normalization volume. 

The atomic transition dipole operator is 



iV 



£ = Ai£(B„+B+). 



(7) 



The matter-field coupling is given formally by the electric 
dipole interaction Hi = ft ■ E. In the rotating wave 
approximation and for linear polarization, we get 



Hr 



qA,n 



T,\I^J7^- e ^ 



2e V 



x [a^Bl e^ na - a\ x B n e-^ na Y (8) 

In the following we treat a single electronic excitation in 
the atomic chain. We start in treating the most symmet- 
ric collective state and then the general collective state. 



A. Symmetric Collective Excitation 

We consider a single excitation in the system with the 
symmetric collective state 



"-^S 



\gi,--- ,e h --- ,g N ). 



(9) 



This state is an eigenstate of the Hamiltonian in the limit 
of qA& > 1 with J /T a < 1, where the atoms are almost 
independent. The other limit of qacl < 1 treated by us 
in other work flll4l8| . 

We calculate the damping rate of such collective state 
through the emission of a photon into free space and the 
damping into the final ground state 



l/> = \gi, 



,9n}- 



(10) 



We apply the Fermi golden rule to calculate the collective 
symmetric state damping rate 



2tt 



^2\(f\Hi\i)\ 2 d(E A -E ph ), (11) 



qA 



which in the present case reads 



Y^ ircq 2 



qA 



e VN 



J2e-^ 



Tl=l 



5(E A - Eph). 
(12) 



The summation over the photon polarization yields 



v^/ n2 2 (q-m) 

2^(/x-e qA ) =/« — 

\ 1 



(13) 



The summation over q can be converted into the integral 



E 

q 
We use 



V 



(2tt) 3 



dO sin 9 I q dq. 

il Jo 



q = g(sin 9 cos (f>, sin 9 sin </), cos ( 
and the transition dipole is taken to be 
/x = /i(sin</2, 0, cost/?). 



(14) 



(15) 



(16) 



The integration over <f>, and the change of the variable 
U = qA® cos9, gives 



r, = - 


,,2 2 f+qAa 

^ qA dy 
iireohaN J- qAa 


iV 
n=l 




X 


v 2 

(1 + cos 2 ip) - . y ., (3 cos 2 tp - 1) 

(qAa) 2 


sing the relation 




JV 


ny 


2 N 

= N + £ S 

n<m— 1 


!cos[(n — 


nM, 



(17) 



(18) 



we reach, after the integration over y, the result 

r, - r A 

where 



JV 



1 + A ^J F[« A o(n-m)] 

n< ?7i—l 



(19) 



n*) - N^(i-cosV) 



cos a; sin x 



l-3cos 2 </>) 



(20) 



In figure (5) we plot the function F(x), for two different 
polarization directions. Using the previous numbers, for 
p = 0° we get F(0.5) = 0.9752, and for p = 90° we get 
0.9507, which justifies the use of F = 1 for optical lattice 
ultracold atoms in our previous works |lll - ll8J . 



LL 




FIG. 5: The function F(x) vs. x. The full line is for p = 0°, 
and the dashed line for <p = 90°. 

For different number of atoms we get 

r.(i) = r A , 

r,(2) = r A {l + F(g A a)}, 



T s (3) r A <j 1 + - [2F( 9A a) + F(2q A a) 



(21) 



The symmetric damping rate can be written in the form 

C N-l 



T s {N)=T A l 1 + 2^ 



(N-n) 

~N 



F(q A an) 



(22) 



In figure (6) we plot T/T A as a function of q A a for TV = 5, 
and for the polarizations ip = 0° and ip = 90°. 

Lets consider F{q A an) to represent a bond between 
two atoms that separated by a distance (an), then in the 
above summation the function F(q A an) is multiplied by 
the number of bonds of this length which is (TV — n). In 
the limit of q A a «1 we get F ~ 1, and then r s (TV) w 
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FIG. 6: The symmetric scaled damping T/Ta vs. qAa, for 
TV = 5. The full line is for ip — 0° , and the dashed line for 
1/2 = 90°. 



iVr^. In the limit of q A a ^> 1 we get F ~ 0, and then 
T S (N) « r^4 with oscillations. 

Now we emphasize the dependence of the symmetric 
state damping rate as a function of the number of atoms 
N. We plot the scaled damping rate T S /T A as a function 
of N for different values of q A a. In figures (7 — 9) we 
plot for q A a — 0.001, q A a = 0.1, and q A a = 1, in the 
two cases of <p = 0° and p = 90°. The damping rate of 
the symmetric state grows linearly with the number of 
atoms for small N, and approach a finite value for large 
N. For q A a > 1 the damping rate approach the finite 
value faster than for q A a <C 1. In figure (10) we plot 
T S /T A as a function of <f> for N — 100 at q A a = 0.1. 
Significant difference appears between the damping rates 
for <p = 0° and <p = 90°. 
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FIG. 7: The symmetric scaled damping rate T/Fa vs. TV, for 
qAa = 0.001. The full line is for ip — 0°, and the dashed line 
for 93 = 90°. 
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FIG. 8: The symmetric scaled damping rate F/Fa vs. N, for 
qAa — 0.1. The full line is for ip — 0°, and the dashed line for 
V? = 90°. 
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FIG. 10: The symmetric scaled damping rate F/Fa vs. ip, for 
AT = 100 at q A a = 0.1. 
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FIG. 9: The symmetric scaled damping rate F/Fa vs. N, for 
qAa = 1. The full line is for tp — 0°, and the dashed line for 
V? = 90°. 



B. General Collective Excitation 

Here we consider the case of a single excitation but for 
a general collective state, which is given by 

\i) = -Jjr^ C *\9ir-- ,e;,--- ,9n), (23) 

where C« = ±1, and J2i Cf = N. Equation (fI7|l reads 



2 2 

VI A 

8ir€ohaN 



dy 



N 



E c - 



-iny 



(1 + cos 2 p) - — !L_ (3 cos 2 tp-l) 
{qAa) 2 



(24) 



We use 

N 



£c„6 



n=l 



AT 



= ^+ ]T C„C m e-'( n - m ^, (25) 

n^m— 1 



where C n Cm — ±1. After integration over y, we get 
r(JV) = T A J 1 + 1 ^ C„(7 m F[g A a(n - m)] 1 



n<m— 1 



(26) 

Here we present the results for two examples. For N = 
2 we have the symmetric state 



,., |ei,g2) + |gi,e 2 ) 

|z)s = 7=2 ^ 

with the damping rate 

r s =T A {l + F(q A a)}, 
and the antisymmetric state 

,.\ |ei,g 2 ) - |.gi,e 2 ) 

|i>Q = 71 

with the damping rate 

T a = T A {l-F{q A a)}. 



(27) 



(28) 



(29) 



(30) 



The results for N — 2 agree with the known results [22j. 
In figure (11) we plot T/T A for the symmetric state of 
N = 2 as a function of q A a, for the polarizations ip = 0° 
and </3 = 90°. In figure (12) the plot is for the antisym- 
metric state. 

For N = 3, for the symmetric state 



l'<>. 



\ei,92,93) + |gi,e 2 ,g 3 ) + |gi,ff2,e 3 ) 



(31) 
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FIG. 11: The scaled damping Y /Ta vs. qAa, for N — 2 with 
the symmetric state. The full line is for p = 0°, and the 
dashed line for p = 90°. 
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FIG. 12: The scaled damping r/r A vs. g A a, for iV = 2 with 
the antisymmetric state. The full line is for <p — 0°, and the 
dashed line for y> = 90°. 



we get 



r> = r A { 1 + - [2F(q A a) + F(2q A a)] 



For the antisymmetric state 

,.\ |ei,52,33) - |.gi,e 2 ,53> + |5i,32,e 3 ) 



>/3 



we get 



r a = T A \ 1 - - [2F(q A a) - F(2q A a) 



(32) 



(33) 



(34) 



In the limit of q A a <C 1 we have F ~ 1, then for the sym- 
metric state we get T ss 3r^4, and for the antisymmetric 
one we get T rj T^/3. In the limit of q A a 3> 1 we have 



F ~ 0, then for the symmetric and antisymmetric states 

we get r « r^. 

In hgurcs (13) we plot T/T A for the symmetric state of 
N = 3 as a function of g^a, for the polarizations tp = 0° 
and (/? = 90°. In figure (14) the plot is for the antisym- 
metric state. 
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FIG. 13: The scaled damping T /Ta vs. qAa, of the symmetric 
state for N — 3. The full line is for tp — 0°, and the dashed 
line for p — 90°. 
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FIG. 14: The scaled damping V /Ta vs. g^a, of the antisym- 
metric state for N — 3. The full line is for <£ = 0°, and the 
dashed line for tp = 90°. 



IV. COLLECTIVE EXCITATION EMISSION 
PATTERN 

Here we calculate the emission pattern of a collec- 
tive state in a chain of N atoms separated by a dis- 
tance a. The transition dipolc of each atom is \X = 
fi(smtp, 0, cos tp), at positions R n = (0, 0,R n ). For sim- 
plicity the observation point is taken to be at r = (x, 0, 0), 



as seen in figure (15). We concentrate here in the limit of 
qA<i > 1 where the atoms can be treated independently. 
The other limit of qao, < 1 investigated by us in previous 
work [18[ . The positive electric field operator of the atom 
(n), in the far zone field where x~^> \a, is given by [24j 



K +> (r,t) 



fiq A sin <pn 

47re |r - Rn 



B t- 



|r-RJ 



e„, (35) 



where <p n is the angle between \x and r — R„, and the 
unit vector e„ is defined by 



e„ = n„ x y, n n = 



r-R„ 



(36) 



|r-Rn| 
We have 

r - R„ = (x, 0, -R n ), |r - R„| 2 = x 2 + B? n , (37) 
and 

4> n = n - Lp - a„, tana„ = x/R n , (38) 

with 

_ (x,0,-Rn) „ _ (-R»,0,a:) 



V* 2 + ^ 



V^+lf 




FIG. 15: The observation point at r along the x axis, and the 
lattice is along the z axis. The angle between the transition 
dipole fi and r — R„ is <p n . The electric field direction off 
atom (n) is e„. 

For the atomic transition operators we use the expec- 
tation values 

(Bi(t-U)) = (S l (0))e-^(*-*' ) e- r ^( t ~*-)/ 2 , 
{B\{t-U)B t {t-U)) = {B\{Q)Bm)e- TA{t - U \ (40) 

and 

(Btit-UWt-ts)) = (BHo^me-^-^^ 



x e 



-iOJA{t%— tj) 



(41) 



where ij = |r — Rj|/c. In the limit qAa > 1 the single 
excitation collective states decay with the single excited 
atom damping rate Ta- 



The total electric field at the observation point is 

E(+)(r,t)=^E| +) (r,i), (42) 

i 

and the intensity is 

/(r,t) = i eo c(E(-)(r,t).E( + )(r,i)). (43) 

Explicitly we can write 

J(r,t) =£!*(!•,*) + £G y (r,i), (44) 

i ijtj 



where the «-th intensity is 



1 



i\( _ )f,. *-\ *?.(+)/ 



/ i (r,t) = -c c(Er(r J t)EP(r,t)) J (45) 



and the correlation function is 



Gy(r,t) = -e c<E^(r,t).E5 + >(r,i)>. (46) 



We get 



J '(-"- 3^CT (B| "" B ' ( ° ))e - r """' 



(47) 



and 



^ l ' j " 327r2 £oC 3 P 



x (Bt(0)B,-(0)) 



e ~ l v * 3 " J e 
sin (A,- sin < 



^^^^Ir-R.Mr-R,!^- 11 ^' 



(48) 



Two-Atoms Chain 



We present the results for the simple case of two atoms. 
One atom is located at the origin Ri = (0, 0, 0), and the 
second at R2 = (0,0, a). The observation point is at 
r = (x, 0,0), where r — Ri = (x, 0,0), and r — R 2 = 
(x, 0, —a), with |r — Ri| = x, and |r — R 2 | = \Jx 2 + a 2 . 
We have the angles <j)\ = § — ip, and <A 2 = tt — <p — a, 
where tana = x/a. We get the times t\ = x/c, and t<i = 
\Jx 2 + a 2 /c. Also we have the unit vectors ei = (0, 0, 1), 
and e 2 = -fett, then ni = (1,0,0), and n 2 = %4r4 
hence (ni«n 2 ) = , ,f , ,, . We obtain 



/i(r,i) 



V^+a^ 



mM "^ (B!(0)B l( 0)) e-^-t), 



327r 2 eoc 3 x 2 



^ = ^3^(^(0)) 
327r z eoc :S or + a 2 

x e 



(49) 



and 



OMr,t) = 4^ ^^p- <Bt(0)B 9 (0)> 

- TA r{-^^)\ e ~ luJA \ — s — ). 

„M sin^sin^ (B t (0)jBl(0)) 



x e 



G 2 i(r,t) 



x e 



327r 2 eoc 3 a; 2 + a 2 



(50) 



Now we consider the two initial states of symmetric and 
antisymmetric collective states. 
For the symmetric collective state 



\i) 



|ei,g 2 ) + |gi,e 2 ) 
V2 



(51) 



we have 



(i? 1 t (0)i? 1 (0)) = (Bt(0)B 2 (0)) 
= (Bt(0)B 2 (0)) = (B^O^O)) = 1 



(52) 



then we get 
r/ ,x Jo (a;) i . 2 



ar + or 



-r\> t- 



yV^ 



a; sm q>\ sm 02 



-r A t- 



x e 



2 cos 



UlA 



VrT? 



(53) 



where we defined the intensity 
Io(x) 



2 4 



167r 2 eoc 3 a; 






(54) 



In figures (16 — 18) we plot the relative intensity 
I(r,t)/Io{x) as a function of a for the angles ip = 0°, 
if = 45° and if — 90°, at the observation point x = 10 6 A 
at the moment i = 2x/c. We use 5a = 1 eV, /i = 1 eA 
and T A = 10 8 Hz. 

For small a the relative intensity is maximum for the 
polarization angle if = 0° and decreases for larger angles. 
It is half for if — 45°, and becomes zero for ip = 90°. The 
maximum of the relative intensity moves into larger a 
with increasing the angle ip. The relative intensity oscil- 
lates in changing a and tend to a finite value for large 
a. Interesting case is for if = 90°, where the intensity 
is zero for small a and increases with increasing a till it 
reach a maximum at a = 10 6 A (for the given numbers), 
and decreases back towards a finite value for larger a. 



In the limit of x >• a, where V ' x 2 + a 2 ~ x + |-, and 
as (f> w -| — </?, we can write 

'('•«> «s^?«- r ' ( " ,) ~ > i * 

<j i + e r -45k +2cos Li— J e TA ^ M55) 
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FIG. 16: The symmetric state scaled intensity I(r,t)/Io(x) 
vs. a, for yj = 0° at x — 10 6 A and t = 2a;/c. 
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FIG. 17: The symmetric state scaled intensity I(r,i)/Io(x) 

vs. a, for ip = 45° at x = 10 6 A and t = 2x/c. 



For the antisymmetric collective state 
|ei,32> ~ |ffi,e 2 ) 



\i) 



V2 



(56) 



we have 



(Bj(0)B!(0)) = (51(0)52(0)) 
(51(0)52(0)) = (B|(0)B 1 (0)) 
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FIG. 18: The symmetric or antisymmetric state scaled in- 
tensity I(r,t)/Io(x) vs. a, for <p = 90° at x — 10 6 A and 
t = 2xj c. 



then we can write 

j i *\ 7 °( x ) / • > 
I(r, t) = — ; — < sin 



„2„- 2 



4 
-r A l t- 



x e 



-rx(t-f) , g sin <h 

1 loo 

ar + a 2 



x 2 sin <j)i sin 02 



2 cos 



x e 



x 2 + a 2 



LOA 



x — y/x 2 + a 2 



(58) 



In figures (19 — 20) we plot the relative intensity 
I(r, t)/Io(x) as a function of a for the angles <p = 0° and 
99 = 45°. The case of cp = 90° is the same as in figure 
(23). As before, the observation point is at x = 10 6 A at 
the moment t = 2x/c, with the other previous numbers. 
The results are similar to the symmetric ones except from 
the case of small a where the relative intensity tends to 
zero as expected. 

In the limit of x 3> a, as (/> « ^ — ip, we can write 



I(r,t) 



2 4 

(J-Ua 



647r 2 eoc 3 x 
1 



., . r -(*-t) C os 2 ^ 



e r ^fe 



2cos( w A ) e rA fe 1.(59) 



V. SUMMARY 

In the present paper we investigated optical properties 
of a one dimensional atomic chain, in which the lattice 
constant can range from a few angstroms up to thou- 
sands of angstroms. Namely, the lattice constant can 
change from being smaller than the atomic transition 
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FIG. 19: The antisymmetric state scaled intensity 
I(r,t)/I (x) vs. a, for ip = 0° at x = 10 6 2 and t = 2x/c. 
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FIG. 20: The antisymmetric state scaled intensity 
J(r, t)/I (x) vs. a, for ^ = 45° at x = 10 6 A and t = 2x/c. 



wavelength up to much larger one. In the limit of lattice 
constant smaller than the atomic transition wavelength 
the electrostatic interactions are applicable, which found 
useful for most of the typical experiments on optical lat- 
tice ultracold atoms. In our previous work we limited the 
discussion to electrostatic interactions, where we consid- 
ered only resonance dipole-dipole interactions, and which 
is justified in the present work. For small lattice constant 
the electrostatic interactions are responsible for the for- 
mation of excitons, where we did extensive study in this 
regime with emphasize on the exciton life times. For large 
lattice constant the inclusion of radiative corrections are 
necessary, which is the main issue in the present paper. 
For large lattice constant the radiative corrections are 
included, and in this regime we found that the coupling 
parameter for the energy transfer among even the near- 
est neighbor atom sites is smaller than a single excited 
atom damping rate. Hence, the energy transfer is not 



10 



favorable, and we treated the atoms as independently 
setting on the lattice sites. Then, we calculated the 
damping rates of different collective electronic excitations 
in including the radiative corrections by considering the 
effect of the existence of all the other atom sites, de- 
spite their large distances from the excited atom. Big 
attention we gave for the most symmetric state, where 
we emphasized the dependence of its damping rate on 
the number of atoms for different lattice constant. We 
found the symmetric damping rate to behave linearly at 
small atom numbers and saturate at large numbers. The 
damping rate of symmetric and antisymmetric collective 
states tend to that of a single excited atom with oscilla- 
tions due to the radiative effect through the exchange of 
virtual photons. The differences between damping rates 
of collective states appear for small lattice constant, in 
which the symmetric states have superradiant damping 
rate, that is N time the single excited atom rate. Here, 
part of the antisymmetric states become dark with zero 
damping rate, and other part is metastable with a frac- 
tion of the single excited atom damping rate. Moreover 
we calculated the emission pattern off a chain of atoms 



with a large lattice constant in which the atoms can be 
considered independently. The emission intensities off 
two atoms with symmetric and antisymmetric states are 
presented as a function of the interatomic distance. 

The results of the present paper are illustrated in terms 
of optical lattice ultracold atoms, but they are general 
and can be adopted for any chain of optically active ma- 
terial. For example, chains of semiconductor quantum 
dots fit exactly in the regime of large lattice constant, 
where radiative corrections are unavoidable, and the life 
times of their collective states can be treated according to 
the present paper. Other system that exploits the regime 
of the present paper is a lattice of large organic molecules 
sitting on a matrix with a given large lattice constant, 
the collective damping rate and emission pattern are ex- 
pected to behave according to our present results. 
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